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Abstract 



This is a summary of the contributions on the next-to-leading order corrections to the BFKL 
equation which were presented to the 'Small- a: and Diffraction' working group at the 1998 Durham 
' Workshop on HERA Physics. 



— 1 < i 

^ ■ The original BFKL equation |l| allowed the computation of elastic scattering amplitudes in the 



Regge limit to leading logarithmic (LO) accuracy, i.e. all terms ~ (a s lns) n are summed. The key 
element is the gluon Green function: 

F(s ' fc)= /^2^(^) u - a sX (j) (fcf) ■ (1) 
The scale so is arbitrary, a s = N c a s /-K is fixed and, at leading order, 

x(7)=Xo(7) = -27B-^(7)-V>(l-7)- (2) 

The s — > oo behaviour of the amplitude is governed by the right-most w-plane singularity which occurs 
at 7 = 7 where 

w = a s x(7)- (3) 

The conformal invariance of the LO kernel allows us to write equation (||). The NLO corrections Q 
break the conformal invariance by leading to the running of the QCD coupling. However, it is possible 



to study the conformally invariant part of the NLO equation. In this case, equation (||) is still relevant 
but now the kernel is 

X(l) = Xo(t) + OL a Xl(i). (4) 

The NLO kernel depends upon the scale so i n a straightforward way (sq = k\k2 is usually chosen). 
The NLO corrections contained in xi are huge for typical values of a s which is a cause for serious 
concern since we are attempting a perturbative expansion ^, |[. For example, at 7 = 1/2 

X(7) = ^o(l - 6.5 a s ) 

where ujq = a s xo(l/2) = 41n2 a s is the leading eigenvalue of the LO BFKL kernel. 

Marcello Ciafaloni and Gavin Salam pointed out that the problem is alleviated by identifying 
and resumming a large part of the NLO corrections: that with double logarithms in the transverse 
momenta (corresponding to the I/7 3 and 1/(1 — 7) 3 parts of the NLO kernel). After ensuring the 
consistency of these double logarithms to all orders, the perturbation series is much more convergent. 

The double logarithms are closely associated with the choice of scale sq. Schematically one can 
see that they arise through a change of scale sq from k\ to k\k2- In the situation with k\ S> &2, the 
quantity which exponentiates is roughly 

~ as ln ln (|) = ^ ln (kk) 111 (|) - r s ln2 (I) ■ 

If one writes one's expansion in terms of scale sq = k\k<i (i.e. the RHS) the exponentiation leads 
not just to a large NLO term, but also to a whole series of NNLO, NNNLO, . . . terms with double 
transverse logarithms, which can be large. These should therefore be resummed. At first sight the 
simplest way to do so seems be to rewrite one's expansion in terms of so = fef (LHS). The problem is 
that such a scale choice spoils the convergence (again by double logarithms) of the kernel in the case 
when k\ <C k%, where the 'natural' (DGLAP) scale choice is instead Sq = k%. 

A way of treating both regions on an equal footing while correctly resumming the double logarithms 
was proposed in || and is related to an approach suggested by the Lund group ||. Figure 1 shows that 
such a resummation greatly improves the convergence of the kernel. The various schemes correspond 
to different treatments of less singular terms (in particular those related to collinear branching and 
the running coupling), which remain to be accounted for in such a way as to ensure consistency with 
the renormalisation group (see below, and also 

The phenomenological importance of matching to the known DGLAP limit and of including kine- 
matical and other 'finite z 7 effects, e.g. coherence, have been and will continue to be explored by a 
number of groups. For example, see the contribution of Kharrazhia and Stasto in these proceedings 
||. Lynn Orr showed us the crucial importance of imposing the correct kinematics. She focussed on 
the azimuthal decorrelation in dijet production at the Tevatron where the data were compared with 
results from a BFKL Monte Carlo with the correct kinematics imposed |J. 

So far we ignored the fact that the QCD coupling runs. At NLO we need to account for the 
breaking of conformal invariance due to running coupling effects. The eigenfunctions of the kernel 
are no longer simple powers of the momenta and a more sophisticated approach is required. Since 
the BFKL formalism integrates over all momenta, the running of the coupling invariably leads to 
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Figure 1: The result of resummation on \ an d its second derivative; shown for rij = 0. A sym- 
metric version of the pure NLO kernel has been used to redefine k 2 ^'^ as y / a s (k 2 )k 2 (' y ~ 1 ^ in the 
eigenfunctions. 
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divergences associated with the failure of QCD perturbation theory in the infrared. To regularise the 
theory, one might take a phenomenological approach and introduce by hand some regulator, e.g. a 
gluon mass to freeze the QCD coupling at some low scale. This has been done before, in the case of 
the LO BFKL equation |1Q[| . However, in single scale problems like deep inelastic scattering at small 
x one does not need to speculate. The collinear factorisation of deep inelastic cross-sections allows 
the ill behaved infrared dynamics to be factorised into parton distribution functions and it is their 
scale dependence that is predicted by renormalisation group (RG) (i.e. DGLAP) equations of QCD. 
Ciafaloni stressed the importance of working within a RG consistent framework. 

A systematic approach to small- x structure functions was presented by Robert Thorne. A few 
years ago Catani showed how to write the DGLAP equations in a manifestly factorisation scheme 
independent way |Q: 

UFl (5) 



SlnQ 2 

where i,j = 2, L and is the physical anomalous dimension matrix. The NLO corrections to the 
BFKL equation can be used to determine the scale dependence of the structure functions i*2 and 
In particular, the physical anomalous dimension matrix is computed including the complete two- loop 
DGLAP splitting functions supplemented with the leading twist information coming from the NLO 
BFKL equation; corrections to this are strictly 0(a s (fi)). The only ambiguity in this procedure is in 
the choice of the renormalisation scale /i. Thorne pointed out that the 'natural' choice /i = Q has 
disasterous consequences, as can be seen in the upper of Fig.2: the dotted line is the evolution of the 
longitudinal structure function at low x with /x = Q at NLO, as compared to the LO evolution shown 
by the solid line (in both cases Fl(x,Q 2 ) is parameterized as (0.1/x) -0 ' 3 at ln(Q 2 /A 2 ) = 8). Thorne 
adopted the BLM scale fixing procedure |I2|] which, for some observable 



R(Q) = Axx a (ji) - (3 a s (fi) 2 [Aln(Q/fi) + B] + a s (p) 2 C + • • • , 
states that the appropriate scale is determined by ensuring that 

A]n(Q/fx)+B = 0. (6) 
In this case A and B depend upon a s and x, and the scale choice corresponds to 



A)(hi(Q 2 /A 2 ) + 3.6a s (Q 2 )(ln 1/x) 1 / 2 ) 

at small x and for all physical anomalous dimensions (and parton anomalous dimensions in sensible 
schemes). This scale choice has the effect of absorbing the large higher order corrections into the 
running of the coupling, i.e. subsequent corrections are small and support the use of perturbation 
theory, as shown in the lower of Fig.2. It also has the peculiar effect of reducing the coupling at low x. 
The physical origin of such strange behaviour was the cause for some discussion, though it now seems 
to be due to infrared diffusion in the BFKL equation only affecting the inputs for structure functions, 
while the evolution is sensitive to ultraviolet diffusion, and hence scales larger than Q 2 as x decreases 
|0|. Thorne also fixed the BLM scale for T22 (~ Pqq) at high x, and showed us the result of fitting 
all the DIS data fitted by the MRST group using lowest order evolution with the scale fixing. The 

1 Actually we don't have all the information yet since the coupling of the gluons to the photon via a quark loop has 
not yet been computed at NLO. 
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Evolution of F L , ln(Q 2 /A 2 )=8, F L =(x/0.1)" 3 
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LO+NLO (exponentiated) 




Figure 2: Evolution of F L at LO and NLO. 



resulting fit produces a remarkably small chi squared, as can be seen in Table 1 which compares the 
latest MRST fit with the lowest order BLM fit (LO(x)). It also leads to F L (x, Q 2 ) following the shape 
of F2(x,Q 2 ) at low Q 2 , unlike the conventional approaches. 
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Data set 


No. of 


LO(x) 


MRST 




data pts 






HI ep 


221 


149 


164 


ZEUS ep 


204 


246 


270 


BCDMS up 


174 


241 


249 


NMC HP 


130 


118 


141 


NMC fid 


130 


81 


101 


SLAC ep 


70 


87 


119 


E665 HP 


53 


59 


58 


E665 nd 


53 


61 


61 


CCFR 


66 


57 


93 


CCFR F% N 


66 


65 


68 


NMC ra/p 


163 


176 


187 


total 


1330 


1339 


1511 



Table 1: Comparison of x 2 for fits using the full leading order (including ln(l/x) terms) renormalization 
scheme consistent expression with BLM scale, and the two-loop MRST fit. 
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